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Abstract 

We describe a method of obtaining closed-form complete solutions of cer- 
tain second-order linear partial differential equations with more than two in- 
dependent variables. This method generalizes the classical method of Laplace 
transformations of second-order hyperbolic equations in the plane and is based 
on an idea given by Ulisse Dini in 1902. 



1 Introduction 

Factorization of linear partial differential operators (LPDOs) is often used in modern 
algorithms for solution of the corresponding differential equations. In the last 10 
years a number of new modifications and generalizations of classical algorithms for 
factorization of LPDOs were given (see e.g. [21 [TOl HU H31 HH HZl HH [20]). Such 
results have also close links with the theory of explicitly integrable nonlinear partial 
differential equations, cf. [Tl [L8l 122] . 

As one can see from simple examples (cf. pj5] and Section [2] below) a "naive" 
definition of factorization of a given LPDO L as its representation as a composition 
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L = L\ o L2 of lower-order operators does not enjoy good properties and in general 
is not related to existence of a complete closed-form solution. 

On the other hand, for second-order hyperbolic linear equations in the plane we 
have a well established and deep theory of "generalized" factorization. This theory is 
known since the end of the XVIII century under the name of Laplace cascade method 
or Laplace transformations. As proved in [19], existence of a complete solution of 
a given second-order hyperbolic equation in the plane in explicit form is equivalent 
to some "generalized factorizability" of the corresponding operator which in turn is 
equivalent to fmiteness of the chain of Laplace transformations ending in a "naively" 
factorizable operator. We give a short account of this method in Section [2j 

There were some attempts to generalize Laplace transformations for higher- 
order operators or larger number of independent variables, both in the classical 
time [121 dSl EES] and in the last decade [21 [20] ■ A general definition of generalized 
factorization comprising all known practical methods was given in [TH]. Unfortu- 
nately the theoretical considerations of [19] did not provide any algorithmic way of 
establishing generalized factorizability of a given LPDO. 

In this paper we move a bit further, extending the algorithmic methods for 
generalized factorization to the case of second-order operators in the space of more 
than two independent variables following an approach proposed by Ulisse Dini in 
1902 HE]. 

The paper is organized as follows. In the next Section we give an exposition of 
the classical theory of Laplace transformations. In Section [3] we work out an example 
which demonstrates the idea of Dini transformations and in Section H] we prove a new 
general result showing that this idea gives a practical method applicable to arbitrary 
hyperbolic second-order linear equation in the three-dimensional space provided the 
principal symbol of the operator factors. The last Section is devoted to discussion of 
algorithmic problems encountered in the theory of Laplace and Dini transformations, 
their relations to the given in [TJ5] theoretical basis. Some conjectures on possibility 
of generalized factorization and existence of complete solutions in closed form are 
given. 

2 Laplace and generalized Laplace transformations 

The cascade method of Laplace (also called the method of Laplace transformations) 
is until the present date the most general method of finding closed-form complete 
solutions of hyperbolic second-order linear partial differential equations with two in- 
dependent variables. Here we briefly sketch this classical theory in a form convenient 
for our purpose. The complete account may be found for example in [SJ [9]. 

Let a general second-order linear operator L with two independent variables x 
and y be given: 

2 

L = ^PiDiD 2 - 1 + ai(x, y)D x + a 2 {x, y)D y + c{x, y), (1) 

i=0 
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Pi = Vi{ x iV)i D x = d/dx, D y = d/dy. Hereafter we will always suppose that the 
operator (jTJ) is strictly hyperbolic, i.e. the characteristic equation X 2 po — ^Pi+P2 = 
for the principal symbol of (jTJ) has two distinct real roots Xx(x,y), \ 2 (x,y), so we 
can introduce two first-order characteristic operators Xi = m i (x,y)D x + n i (x 1 y)D y) 
i = 1,2, mi/rii = Aj (Xi are defined up to rescaling Xi — > "fi(x, y)Xi). 

The corresponding equation Lu = now may be rewritten in one of two charac- 
teristic forms: 

(XiX 2 + a x Xi + a 2 X 2 + a 3 )u = 
(XjCx + oixXx + a 2 X 2 + a 3 )u = 0, 

where y). Since the operators Xi do not necessarily commute we have to 

take into consideration in (j2J) and everywhere below the commutation law 

[X u X 2 ] = X X X 2 - X 2 X X = P(x, y)Xx + Q(x, y)X 2 . (3) 

Using the Laplace invariants of the operator (J2J) 

h = X 1 (a 1 ) + a x a 2 — a 3 , k = X 2 (a 2 ) +a{a 2 — a 3 , 

we represent the original operator L in two possible partially factorized forms: 

L — (Xi + a 2 )(X 2 + ai )-h= (X 2 + a^Xx + a 2 ) - k. (4) 

From these forms we see that the equation Lu = is equivalent to any of the 
following two first-order systems 

/CM \ X 2 u = -a x u + v , , ( Xiu = -a 2 u + w, , . 



Xxv = hu — a 2 v. i X 2 w = ku — a\w. 

If at least one of the Laplace invariants h or k vanishes identically, then the operator 
L factors (in the "naive" way) into composition of two first-order operators and the 
corresponding system in (jSJ) becomes triangular. So the problem of integration of 
the original second-order equation is reduced to a much easier problem of integration 
of linear first-order equations. The latter problem is essentially reducible to finding 
the complete solution of a (nonlinear!) ODE, see below Section 

If h 7^ 0, k 7^ 0, one can take one of the systems §J5§ (to fix the notations we 
choose the left system (Si)), express u using the second equation of (Si) 

u = (Xiv + a 2 v)/h (6) 

and substitute this expression into the first equation of (Si) in (jSJ); as the result one 
obtains a X x -transformed equation Ln)V — 0. It has different Laplace invariants (cf. 

DP) 

h [x) = X x (2ax -P)- X 2 (a 2 ) - X X X 2 In h + QX 2 \nh-a 3 + (ax - P)(a 2 - Q) 
= 2h- kXxX 2 \nh + QX 2 \nh + X 2 (Q) - X x (P) + 2PQ, 

k(x) = h. 
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If h(i) = 0, we can solve this new equation in quadratures and using the same 
differential substitution ([HD we obtain the complete solution of the original equation 
Lu = 0. 

If again hn) ^ 0, apply this Xx-transformation several times, obtaining a se- 
quence of second-order operators L( 2 ), 1/(3), ... of the form (j2J). If on any step we 
get h( k ) = 0, we solve the corresponding equation Lqauha = in quadratures and, 
using the differential substitutions (JBj), obtain the complete solution of the original 
equation. Alternatively one may perform X 2 -transformations: rewrite the original 
equation in the form of the right system (S 2 ) in ([5]) and using the substitution 
u = (X 2 u> + a\w)/k obtain the equation Lr-nw = with Laplace invariants 

h-D = k , (7) 

fc ( _i) = 2k - h - X 2 XilnA; - PXilnA; + X 2 {Q) - X X (P) + 2PQ. 

In fact this X 2 -transformation is a reverse of the Xi-transformation up to a gauge 
transformation (see [1]). So we have (infinite in general) chain of second-order 
operators 

^2 r X2 f X2 f X\ i X\ i X\ ro \ 
• ^(-2) L (-l) <— L —> L (l) —> L (2) — ► • • • (oj 

As one may prove (see e.g. [2]) if the chain ([S]) is finite in both directions (i.e. we 
have hrm = 0, h^_K) — for some N > 0, K > 0) one may obtain a quadrature-free 
expression of the general solution of the original equation: 



u 



c F + ciF' + ... + c N F {N) + d G + diG' + ... + dx^G^-^ (9) 



with definite Ci(x,y), di(x,y) and F(x), G(y) - - two arbitrary functions of the 
characteristic variables and vice versa: existence of (a priori not complete) solution 
of the form ([9]) with arbitrary functions F, G of characteristic variables implies 
= 0, /i(_ r ) = for some s < N, r < K — 1. So minimal differential complexity 
of the answer (Q (number of terms in it) is equal to the number of steps necessary 
to obtain vanishing Laplace invariants in the chain ([8]) and consequently naively- 
factorable operators. If (IE|) is finite in one direction only, one can still obtain a 
closed-form expression for the complete solution of the original equation; however, 
it will have one of the free functions F or G inside a quadrature expession. More 
details and complete proofs of these statements may be found in [U [9] for the case 
Xi = D x , X 2 = D y , for the general case cf. [91 p. 30] and []]. 

Example 1. As a straightforward computation shows, for the equation u xy — 
u = ^ ne chain ([8]) is symmetric (h^ = /i(_j_!)) and has length n in either 
direction. So the complexity of the answer ([9]) may be very high and depends on 
some arithmetic properties of the coefficients of the operator L; for the equation 
u xy — (j;+ y )2 u = the chain ([8]) will be infinite unless the constant c = n(n + 1) . 

Recently a generalization of this classical method was given in [20]. In is appli- 
cable to strictly hyperbolic linear equations of arbitrary order with two independent 
variables x, y only. 
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In [1] a simple generalization of Laplace transformations formally applicable 
to some second-order operators in the space of arbitrary dimension was proposed. 
Namely, suppose that such an operator L has its principal symbol 

Sym = ^2a ili2 (x)D Xii D Xi2 

which factors (as a formal polynomial in formal commutative variables D x .) into 
product of two first-order factors: Sym = X\X 2 (now Xj = ^\ bij(x)D Xi ) and 
moreover the complete operator L may be written at least in one of the forms given 
in ([2]). This is very restrictive since the two tangent vectors corresponding to the 
first-order operators Xj no longer span the complete tangent space at a generic point 
(xq). ([3]) is also possible only in the case when these two vectors give an integrable 
two-dimensional distribution of the tangent subplanes in the sense of Frobenius, i.e. 
when one can make a change of the independent variables (x) such that Xj become 
parallel to the coordinate plane (x\,x 2 ); thus in fact we have an operator L with 
only D Xl , D X2 in it and we have got no really significant generalization of the Laplace 
method. If one has only §2§ but ([3]) does not hold one can not perform more that 
one step in the Laplace chain (JHJ) and there is no possibility to get an operator with 
a zero Laplace invariant (so naively factorizable and solvable). 

In the next section we demonstrate, following an approach proposed by U. Dini 
in another paper [5], that one can find a better analogue of Laplace transformations 
for the case when the dimension of the underlying space of independent variables 
is greater than two. Another particular special transformation was also proposed 
in [2J, [2T]; it is applicable to systems whose order coincides with the number of 
independent variables. The results of [2J, [21] lie beyond the scope of this paper. 

3 Dini transformation: an example 

Let us take the following equation: 

Lu = (D x D y + xD x D z - D z )u = 0. (10) 

It has three independent derivatives D x , D y , D z , so the Laplace method is not 
applicable. On the other hand its principal symbol splits into product of two first- 
order factors: £i£ 2 + ^£163 = £1(^2 + ^£3)- This is no longer a typical case for 
hyperbolic operators in dimension 3; we will use this special feature introducing 
two characteristic operators X\ = D x , X 2 = D y + xD z . We have again a nontrivial 
commutator [Xi, X 2 \ — D z — X s . The three operators Xi span the complete tangent 
space in every point (x, y, z). Using them one can represent the original second-order 
operator in one of two partially factorized forms: 

L = X 2 X\ — X3 = X\X 2 — 2X3. 
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Let us use the first one and transform the equation into a system of two first-order 
equations: 

Lu = <=► l ±lU = V 2 (11) 
L X 3 u = X 2 v. 

Here comes the difference with the classical case dim = 2: we can not express u as 
we did in (J6]). But we have another obvious possibility instead: cross-differentiating 
the left hand sides of (TIT!) and using the obvious identity [X 1; X 3 ] = [D X ,D Z ] = 
we get X x X 2 v = D x (D y + xD z )v = X 3 v = D z v or = D x (D y + xD z )v — D z v = 

(b x b y + xb x b z ) v = (b y + xb z )b xV = x 2 x lV . 

This is precisely the procedure proposed by Dini in [5]. Since it results now 
in another second-order equation which is "naively" factorizable we easily find its 
complete solution: 

4>(x, xy — z) dx + ip(y, z) 



where and ip are two arbitrary functions of two variables each; they give the 
general solutions of the equations X 2 (j) = 0, Xiip = 0. 
Now we can find u: 



u= (^v dx + (by + xD z )v dz^j + 9(y), 



where an extra free function 9 of one variable appears as a result of integration in 

m 

So we have seen that such Dini transformations (1111) in some cases may produce a 
complete solution in explicit form for a non-trivial three-dimensional equation ffTOl) . 
This explicit solution can be used to solve initial value problems for ffTOl) . 

4 Dini transformation: a general result for dim = 

3, ord = 2 

Dini did not give any general statement on the range of applicability of his trick. In 
this section we investigate this question. Obviously one can make different transfor- 
mations similar to that demonstrated in the previous section, here we concentrate on 
the simplest case of second-order linear equations with three independent variables 
whose principal symbol factors. 

Theorem 1 Let an operator L = J2i+j+k<2 a ijk( x i V-, z )D l x b y D^ has a factorizable 
principal symbol: J2i+j+k=2 a ijk( x > V, z )b x b y D^ = SiS 2 (modulo lower-order opera- 
tors) with (non- commuting) first-order operators Si, S 2 ; S\ ^ \(x,y,z)S 2 . Then in 
the generic case there exist two Dini transformations Lm, £(-i) of L. 

Proof. One can represent L in two possible ways: 

L = SiS 2 + f + a(x, y, z) = S^i + U + a(x, y, z) (12) 
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with some first-order operators T, U. We will consider the first one obtaining a 
transformation of L into an operator Lm of similar form. 

In the generic case the operators Si, S 2 , T span the complete 3-dimensional 
tangent space in a generic point (x, y, z). Precisely this requirement will be assumed 
to hold hereafter; operators L with this property will be called generic. 

Let us fix the coefficients in the expansions of the following commutators: 

[S 2 , f ] = K(x, y, z)St + M(x, y, z)S 2 + N(x, y, z)f. (13) 

[Si, S 2 ] = P(x, y, z)Si + Q(x, y, z)S 2 + R(x, y, z)f. (14) 

First we try to represent the operator in a partially factorized form: L = (Si + 

a) (S 2 + ft) + V + b(x,y,z) with some indefinite a = a(x,y,z), ft = ft(x,y,z) and 
V = f-(3S 1 - aS 2 , b=a-aft- Si(ft). 

Then introducing v = (S 2 + (3)u we get the corresponding first-order system: 

Lu -°^{ (V + b )u=-(S 1 + a)v. (15) 
Next we try to eliminate u by cross-differentiating the left hand sides, which gives 
[(V + b), (S 2 + (3)}u = (S 2 + /3)(Si + a)v + (V + b)v. (16) 

If one wants u to disappear from this new equation one should find out when [(V + 

b) , (S 2 + /3)]u can be transformed into an expression involving only v, i.e. when this 
commutator is a linear combination of just two expressions (S 2 + (3) and (V + b): 

[(V + b), (S 2 + 13)) = fi(x, y, z)(S 2 + (3) + u{x, y, z){V + b). (17) 

This is possible to achieve choosing the free functions a(x, y, z), (3(x, y, z) appropri- 
ately. In fact, expanding the left and right hand sides in (|T7|) in the local basis of 
the initial fixed operators Si, S 2 , T and the zeroth-order operator 1 and collecting 
the coefficients of this expansion, one gets the following system for the unknown 
functions a, /3, /1, v. 

K + (3P-S 2 ((3) = u(3, 
M — S 2 (a) + j3Q = va — ji, 
N + (3R = -v, 

/3S x (/3) - f(J3) + S 2 (a) - /3S 2 (a) - S 2 (Si(^)) = -v{a - a (3 - Si(/3)) - fji/3. 

After elimination of v from its first and third equations we get a first-order non-linear 
partial differential equation for /?: 

S 2 ((3) = P 2 R+(N + P)f3 + K. (18) 

This Riccati-like equation may be transformed into a second-order linear PDE via 
the standard substitution ft = S 2 (j)/j. Taking any non-zero solution ft of this 
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equation and substituting \i = va + 5*2 (a) —(3Q — M (taken from the second equation 
of the system) into the fourth equation of the system we obtain a first-order linear 
partial differential equation for a with the first-order term f3S 2 (a). Any solution 
of this equation will give the necessary value of a. Now we can substitute [(V + 
b), (S 2 + (3)}u = fi(S 2 + P)u + v(V + b)u = fiv — u(Si + a)v into the left hand side 
of (|T6l) obtaining the transformed equation L^v = 0. 

If we would start the same procedure using the second partial factorization in 
(|T2l) we would find the other transformed equation Lr_x\w = 0. □ 

As a rule neither of the obtained new operators £(i), £(-i) factors into a product 
of first-order operators as was the case for the operator L = (D x D y + xD x D z — D z ) 
in the previous section. Then one can repeat the described process due to the fact 
that the principal symbol of the transformed equations is still SiS 2 . Thus we have 
in the case treated in this section an infinite chain of Dini transformations 

• £(-2) ^(-i) *~ L ~~ ^(i) ~^ -^(2) 

If some of the Lu\ is factorizable we can obtain its complete solution (under the 
assumption that one can solve the corresponding first-order equations explicitly) 
and solving the system f|T5l) w.r.t. u step through this chain back (this again requires 
solution of linear first-order equations) finally obtaining the complete solution of the 
original equation Lu = 0. 



5 Open problems 

In the previous three sections we tacitly assumed that the problem of solution of 
first-order linear equations 

Y,bi(x)D Xi +b (x))u = (19) 



or 



bi(x)b Xi + b (x)j u = 0(f) (20) 

can be solved at least for polynomial bi(x). In fact it is well known that even in the 
case dim = 2 the problem of complete solution of (bi(x, y)D x + b 2 (x, y)D y )u = is 
equivalent to finding a nontrivial conservation law for the corresponding nonlinear 
autonomous ODE system or a non-autonomous first-order ODE: 

H = bl(x ' y) ' ^ dy _ b 2 (x,y) 
dy dx h(x,y) 
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(or, equivalently, finding their general solutions). For polynomial bi(x,y) this is 
one of the famous fields of research: study of polynomial vector fields in the plane. 
Recently an essential advance was made in [3J El CTJ; one may hope that a complete 
algorithm may be found. Still the problem of finding complete solutions of (|T9|) in a 
suitable "constructive differential field algorithmically is a challenging problem, as 
well as the problem of finding solutions for the equation ([I8]) . 

Another challenging problem is to establish a connection between the general 
theoretic definition given in [19] and the exposed above practical methods based of 
Laplace and Dini transformations. The known cases suggest the following conjec- 
tures presumably valid for operators of arbitrary order and any number of indepen- 
dent variables: 

• Conjecture 1. If a LPDO is factorizable in the generalized sense of [19], then 
its principal symbol is factorizable as a commutative polynomial in formal 
variables D Xi . 

• Conjecture 2. If a LPDO of order n is solvable (i.e. the corresponding linear 
homogeneous equation has an explicit closed-form solution) then its principal 
symbol splits into product of n linear factors. 

One may also suggest to define the principal symbol of a LPDO using different 
weights for different D Xi ; this would imply for example generalized irreduciblity of 
parabolic operators similar to D 2 X — D y and potentially provide a powerful criterion 
of (un) solvability. 

One should point out that the methods for solution of LPDEs given in the previ- 
ous sections can not be called completely algorithmic: even for the classical case of 
Laplace transformations and the simplest possible characteristic operators X\ = D x , 
X2 = D y we do not have any bound on the number of steps in the chain ([8]). Ex- 
ample 1 given in Section [2] suggests that such bounds or other hypothetic stopping 
criteria would depend on rather fine arithmetic properties of the coefficients of LP- 
DOs. 

A more general theoretic treatment suitable for arbitrary (even under- or over- 
determined systems, cf. [131 [H]) based on the language of abelian categories will be 
exposed in a later publication. 

A link to the theory of Darboux-integrable nonlinear PDEs established in [Tj, 
[T8| [22] in our opinion can be extended to other types on nonlinear PDEs. In this 
connection a generalization of Laplace invariants for higher-dimensional and higher- 
order cases started in [2], [21], [20], [17] would be of extreme importance. 
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